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Abstract 



We investigate the method to combine the techniques of quantum chemisty 
with QED. In our theory, we treat the N-electron system and the Dirac sea on 
an equal footing; we regard both of them as the dynamical degrees of freedom 
of a many-body system. After the introduction of our QED-SCF method, the 
QED-SCF solutions are classified into several classes on the basis of group- 
theoretical operations such as time reversal, parity and 0(3) rotational sym- 
metry. The natural orbitals of general QED-SCF solutions are determined 
by diagonalizing the first order density matrix. Thus, we obtain the possibil- 
ity to treat the system under strong QED effect by the methods of quantum 
chemistry, such as QED-MCSCF and QED-coupled-cluster approaches, 
keywords; QED, QED-SCF, MCSCF, Coupled- cluster theory, Radiation chemistry. 

I. INTRODUCTION 

The Hartree-Fock ( HF ) and multiconfiguration HF ( MCHF ) methods have been well- 
established for nonrelativistic Hamiltonian, and the self-consistent field ( SCF ) equations 
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to obtain the HF and MCHF solutions are also derived in several ways [1]. The most general 
HF solution for the nonrelativistic case is given by the general spin orbitals ( GSO ): The 
two-component spinor [2,3], 



(1) 



(pi = 4>Ja + (piP = 

Very recently we have developed the program package for ab initio HF and density functional 
( DFT ) calculations of molecules by using the GSO [4,5]. It was shown that GSO HF 
and GSO DFT approaches are useful for qualitative understanding of correlation and spin 
correlation effects in polyradicals with orbital degeneracies [4,5]. 

By the way, it is a well-known fact that, in atomic and moleculer physics, the relativistic 
effects become crucial in many cases. Especially in the electronic structure calculation for 
heavy elements, we have to introduce a kind of relativistic treatment. For example, the case 
of neutral heavy elements can be treated satisfactorily by the Dirac-Coulomb-Breit no-sea 
scheme [6]. This method is based on the no-sea (the approximation completely neglects 
the effects of the Dirac sea) Hamiltonian. It has a Dirac-type one-body Hamiltonian, and 
the interaction of electrons is treated via the Coulomb potential and Breit operator [7]. 
This Hamiltonian with several many-body schemes (HF, relativistic many-body perturbation 
theory ( RMBPT ), so on) give good numerical results for neutral heavy elements in atoms 
and molecules. But, especially in the case of highly charged heavy elements, the accurate 
second-order RMBPT indicates that quantum electrodynamical ( QED ) corrections are 
clearly seen in the difference between RMBPT and experiment. Various relativistic schemes 
have already been presented in the electronic structure theory, and they are summarized in 
recent review articles [6,8,9]. 

In the previous paper [10], we constructed the theory of quantum electrodynamical self- 
consistent fields ( QED-SCF ). We derived the time- dependent Hartree-Fock ( TDHF ) 
theory, HF condition and random phase approximation ( RPA ) [11,12]. Because we con- 
structed our theory based on QED, we regard them as an intrinsic treatment of relativistic 
theory in atomic and moleculer physics. In this theory, we treat an N-electron system and 



the Dirac sea on an equal footing, and they interact by exchanging photons. Thus our theory 
can treat the QED effects such as the vacuum polarization effect, which will be observed 
around the nucleus of heavy element ( strong external Coulomb field ) [8,9,13]. To describe 
the QED effects, we have to treat not only N-electron system but also the Dirac sea as the 
dynamical degrees of freedom in a many-body system. Futhermore, because we gain the 
ability to introduce various many-body techniques of atoms and molecules to QED, we can 
discuss the availabilities of post-HF scheme through consideration of the QED-HF stability 
condition, like the GSO model in the nonrelativistic theory. 

This paper is organized as follows. In Sec. II, we give a brief summary of our method to 
combine many-body techniques to QED. In Sec. Ill, we perform group-theoretical character- 
izations of QED-SCF solutions [10] on the basis of time-reversal, parity and 0(3) rotational 
symmetry. Under these group-operations, the first-order density matrix of these solutions 
is classified into several different classes. In Sec. IV, we give some examples to apply our 
method to quantum chemistry. We examine the post QED-SCF methods as in the case of 
post HF methods in quantum chemistry. The MBPT, coupled-cluster ( CC ) and MCSCF 
methods are considered in the QED scheme. To this end, we utilize previous procedures 
for the nonrelativistic case [14]: (1) Use of density matrix, (2) its diagonalization to obtain 
the natural orbitals given by four-component bispinor, and (3) selection of active space for 
MCSCF. Finally in Sec. V, we give discussion about this work. 

II. QED-HF THEORY 

In this section, we give a brief summary of our QED-SCF method which has been inves- 
tigated in previous paper [10]. Our QED-SCF method is based on three factors: First, we 
introduce a QED Hamiltonian written by products of the one-particle operators ( creation- 
annihilation operators ). Second, we introduce an assumption for evaluation of matrix 
elements. Third, we use the Thouless form for the determinantal state [11,15]. By intro- 
ducing these three tools, it becomes possible to combine various many-body techniques with 
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QED. 



A. Hamiltonian 



We start from the following Lagrangian 

C(x) = ~F^(x)F^(x) + *P(x)(iYd^ - mo)^(x) + e$(x)^(x)(A^x) + Af(x)). (2) 

Here, is a ( virtual ) photon field which describes the interaction between fermions, and 
A^ is a classical external field. Next we remove the photon degree of freedom and introduce 
a Coulomb potential U(r) of a nucleus ( or nuclei ) as the external field. Then we obtain 
the following Hamiltonian: 

H = J d 3 xij(x)(-rf ■ V + m + ^°U(r))^(x) 

l -e 2 J d 3 x J d 4 yi>(x)Yip(x)D^(x - y^ivWiv)- (3) 



+ - 



Here each field is quantized. D^ u (x — y) is the full photon propagator. The 7-matrices are 
defined by ( under the standard representation ) 
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where I is the unit matrix, and <jj ( i=l ,2,3 ) is the Pauli matrix: 
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In the above Hamiltonian, we expand each field operator by one-particle state function as 

fa) = £(^ +) (*R + ^ _) (x)SJ), (6) 

i 

k%) = E(^ (+) (^)«i + (7) 

i 

where +(— ) means the electron ( positron ) state, i denotes quantum numbers for one- 
particle states, a|(oj) is electron creation ( annihilation ) operator, S|(Sj) is positron creation 
( annihilation ) operator, and is the four-component bispinor given like 



* (±, = (xf )(1, .X, (±,(2, ,X« ±,,3, ,x! ±),4) ). (8) 
Then we obtain the Hamiltonian written by creation-annihilation: operators 

H = H K + Hi, (9) 
H k = E ^(afo + + - ( 10 ) 
Hi — - V^^ :±± (aja|.a;aj + ajajj&Jaj + a\a,jb k ai + a\b\ajb k 

+a\b]a\<ii + aj^aj.&j' + a\b]b k ai + a\b\b)b k + a^a^ + a^SjSjaj 

+6 i a i 6 fe a; + &]%Mfc + &]4^ + alb]b\bi + b)b k biai + b]b]b k bi), (11) 

where we designate matrix elements as 

7;±± = | rf 3 ^ (±) (x)(- ? 7- V + m + 7ot/(r)Kf } (x), (12) 

(13) 

Here, Ae = je^ — ej| = |e» — e j- j ( the energy difference of one-particle states ). g^ v is 
usual metric tensor and defined as = diag(l, — 1, — 1, — 1). About the interaction, we 
only take into account 0th order bare photon propagator, and choose the Feynman gauge 
for convenience. This Hamiltonian can describes the QED effects such as the vacuum- 
polarization, because which includes the Dirac sea as the dynamical degree of freedom. 
Hereafter we use this Hamiltonian in our theory. 

Next we give the method for evaluation of expectation value. In QED, the Dirac vacuum 
under the presence of the external Coulomb field generates 4-current as an observed effect, 
which is called the vacuum polarization [8,9,13]. It is well-known that, this effect is described 
by the charge conjugation symmetric 4-current given in the next form [8,9,13] 

iTW = -UT, tPM-y^Pix) - £ ft + )(x) 7 ^(x)). (14) 

^n<—m n>—m 
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This expression gives the 4-current induced by the external Coulomb field. In this expression, 
in the external-field-free case two sums cancel, but for the field present case two sums do 
not cancel completely. Hereby, this result should be interpreted as describing the vacuum 
polarization effect which is induced by the external Coulomb field. Therefore, the current 
of N-electron system with adding the vacuum polarization ( here we consider the case that 
the N-electrons occupy the one-particle states up to the Fermi level, as illustrated in Fig. 1 
) is given as 

M*) = -e E $ +) (s)7^ +) (*) 

—m<i<F 

-\{ E $t\xhsi>t\x) - E (is) 

\i<— m i>—m I 

The second line gives the vacuum polarization [8,9,13]. To include this effect in our Hamil- 
tonian, we have to modify the contraction scheme. From the definition of the contraction: 

A'(t)B'(t') = {F\A(t)B(t')\F)6(t - t') - (F\B(t')A(t)\F}6(t' - t). (16) 

Here, A(t) and B(t) are some kind of operators, and |F) is a Fermi sea. We introduce the 
definition for the step function at the same time as 9(t — f)t=t' — Q(f — t)t=v = \- Then we 
obtain the following relation: 

(F\T(A(t)B(t'))\F) t=t , = (F\A-(t)B-(t')\F) t=t , 

= \{F\A{t)B{t)\F) - l -{F\B(t)A(t)\F). (17) 

Let us consider the fact that, it is clear from the Hamiltonian given in Eq. (3), the expecta- 
tion value of it will be written in a functional of the 4-current j^, as same as the discussion in 
Ref. 8. Therefore, we introduce a hypothesis that our Hamiltonian is written by operators 
in the same-time T-products, and when we factorize the vacuum expectation value for an 
operator product with the aid of the Wick theorem, each contraction should be calculated 
by the definition given in (17). 

By using the above method, for example, the expectation value of H is given by 
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(f\h\f) = - £ r«- + -( £ -E)^l + 

i<-m ~m<i<F i>F 

+ E -EX E -E)(^ ++ -^ ++ ) 

L -rn<i<F i>F -m<k<F k>F 

+( E -EX E )(*&T -*&r + ) 

—m<i<F i>F k<—m 

+( E -E)(E)(^ ++ -^ + -) 

-m<k<F k>F i<-m 

+( E )( E )(V£r" " v ikki~~)}- (18) 

i<—m k<—m 

This gives the HF energy in our theory. We argue that we can generalize this method to 
calculate expectation values for arbitrary products of operators. By using the Wick theorem, 
we always factorize the matrix elements of an operator product into the sum of the products 
of contractions, then the contraction is calculated by the definition given above. It is clear 
from the above discussion that only the definition of the contraction for same-time operators 
is modified. 



B. HF solution 



Now, to obtain the relativistic Slater determinant in the Thouless form [11], we will 
obtain a relativistic exponential operator which transforms one to another representations. 
We assume the following relations for the canonical transformation: 
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where S\ = S\ will be satisfied. Here new operators are expanded by a complete set of the 
old representation. We introduce the exponential operator in the following form: 

e lSl = expji E( a «™ a m«n + (3 mn a ] J>i + (3 mn b n d m + 7mn&A)}- ( 23 ) 

mn 
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Here we demand that the matrices formed by the paramerters a mn , (3 mn , (3* mn , ^* mn should be 
Hermitian. Then the exponential operator formally given above is unitary. It is clear that 
the operator (23) will give the relations in (19)~(22). We will write a relativistic Slater 
determinant as 

|$(a,/3,/F,7*)) = e^ QA/r '^|$ >- (24) 
Here the Slater determinant of old representation is expressed by 

l*o> = II«t|0> (25) 

i 

and the Dirac vacuum is defined by Oj|0) = and bi\0) = 0. 

By using the relativistic Slater determinant given above, we can write the expectation 
value of our Hamiltonian into an expanded form, by the same way as nonrelativistic cases. 
The expectation value of our Hamiltonian is given as follows: 

+i(%\[H,S 1 (a,P,P* n *)]\®o) 
+^($ |[[^,5 1 (a,A/5*,7*)],5i(«,A/5*,7*)]|<fo) 
+0(Sf). (26) 

In (26), the first line is the HF total energy, the second line gives the first derivatives with 
respect to the parameters (a, (3, (3*, 7*), and it must be zero in the HF condition. The third 
line corresponds to the second derivatives and they determine the stability of the HF state: 

% ^(* \[[H, Si(a, (3, (3*, 7*)], 5i(a, (3, (3*, 7*)] |$o> > 0, (27) 

and it is equivalent to the stability of collective modes in the RPA. It is obvious fact that the 
operator formalism makes possible to do these discussions in QED. We gave the evidence 
in previous paper [10] that, by using our method, we can derive the TDHF equation, HF 
condition and RPA equation in QED with no inconsistency. Thus we can obtain the most 
stable generalized QED solution at the HF level. 
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III. GROUP-THEORETICAL CLASSIFICATION OF GENERALIZED QED-HF 

SOLUTIONS 



A. Group-theoretical classification 

In this section, we give a brief discussion of group-theoretical classification of the gen- 
eralized QED-HF solutions. We firstly introduce the density matrix like the work of Fuku- 
tome [16]: 

Qaxa = -(i){x)i)(x)) ixi . (28) 

Here ip and ip are usual Dirac field, and they are 4-component bispinors. Thus our density 
matrix is 4x4-matrix, as denoted above. We can expand the 4x4 density matrix into the 
16-dimensional complete set of 7-matrices: 

Q 4X 4 = Q S I + Qli" + + Q A ^Y + Q p i lb . (29) 

In this expansion, we take a convention that S denotes the scalar, V denotes the vector, 
T denotes the 2-rank antisymmetric tensor, A denotes the axial vector and P denotes the 
pseudoscalar. / is the 4x4 unit-matrix, 7^ is usual Dirac gamma matrix, a^ v is defined as 
oiiV = \ [t^jT"]) an d 75 — «7°7 1 7 2 7 3 . Thus we obtain the Lorentz structure in our density 
matrix Q. 

In the relativistic theory, we usually treat the Poincare group ( 4-translation and the 
Lorentz group ), C ( charge conjugation ), P ( parity ), T ( time-reversal ). Under the 
charge conjugation, tp and ip are transformed as 

ij) -> C^ T , $ -> -*fj T C-\ (30) 

Then the density matrix is transformed as 

Q = -(W) -> -C^fC' 1 = CQ T C~ l 

= Q S I - QlY - Q T ^ V + + Q p t l5 . (31) 
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Here C = ry 2 7° is the charge conjugation matrix, and T denotes the transposition of matrix. 
Under the time reversal, 

^(t)-TV(-i), 4>(t)^4,(-t)T, (32) 
together with the rule of taking the complex conjugate about c- numbers, we obtain 
Q(t) = -(mm) -> -T(^-t)i>(-t)YT = TQi-tyr 

= Q s *(-t)I + g y *(-t) 7 ° - QY*(-t)Y + Ql*(-t)* 0i 

-Qj;(-t)av + Q*> (-t)757° " Q^HfrsV " Q P *(-^)^ 5 . (33) 
Here T = ry^ 3 and i, j = 1, 2, 3. Under the spatial inversion, 

^(x) - 7 V(-x), ^(x) - ^(-x) 7 °, (34) 
the density matrix is transformed as 

Q(x) = -Mx$(x)) - -7 (^(-x)^(-x))7 = 7°Q(-x)7° 

= q 5 (-x)/ + g y (-x) 7 ° - QYi-xW - g^(-x)^ 

+Q5(-x)a« - Q^(-x)757° + Qt (-x) 7 57 i " Q P (-x)*75- (35) 

The QED-HF solutions can be group-theoretically classified into several types as in 
the case of nonrelativistic HF solutions [16]. To consider this problem, we determine the 
symmetry-group of a system. In the atomic or molecular systems, the translation invariance 
is broken. In the case of an atom, only 0(3) rotation is remained in the Lorentz group ( 
In the case of a molecule, 0(3) is replaced by the point group.). Under the 0(3) rotational 
symmetry, we expand the Q s , Q p , Qq and Q$ by the scalar spherical harmonics, while we 
expand the Qj , Qf and by the vector spherical harmonics [17]. It is clear from (33) 
and (35), the behavior of each type of the density matrix given in (29) under the spatial 
inversion and time reversal depends not only on the structure of the 7-matrix, but also on 
the angular momentum of the spherical harmonics. Let us consider the case of an atom. We 
treat the group G = 0(3) x P x T. Here, we introduce the subgroup of G as 
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0(3) x PxT, 0(3) x P, 0(3) x T, 

PxT, 0(3), P, T, 1. (36) 

For example, 0(3) x P x T-invariant solution is given as 

Q S * = Q S , Ql* = Q%, others = 0. (37) 

There is no case for 0(3) x P-invariant solution. 0(3) x T-invariant solution is given as 

Q S * = Q S , Qo* = Q%, Qt = Qt others = 0. (38) 

0(3)-invariant solution is given as 

Q S * = Q S , Q%* = Q%, Qo* = Qt Q P * = Q P , others = 0. (39) 

Due to the rotational symmetry, each density matrix can only have an s-wave component, 
in all cases (37)~(39). All solutions given above are for closed shell states. We may have 
magnetic QED-HF solutions for systems under discussions, if time reversal symmetry is 
broken ( For example, the cases of (38) and (39) ). 



B. Nonrelativistic limit 



To consider the nonrelativistic limit, we take the standard representation. The four 
component spinor ip is partitioned as [17] 

(A 



1> 



where <f> is the large component and x is the small component. Then, 
Q s = ^trQ = i(^) = i(0V_ x t x)) 



v 



1 1-1 

Q P = ^tri7 5 Q = -(ifn'ysip) = -HxV + ^x)- 



(40) 



(41) 
(42) 
(43) 
(44) 
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Therfore, at the nonrelativistic limit, the Q s and Qq coincide with each other, while the 
Qq and Q p vanish: Note that, on the other hand, Q s and Q p vanish at the ultrarelativistic 
limit. 

Under the presence of the vectorial density matrices like the Qf , Qf and Q^, the 0(3) 
rotational symmetry will be broken. For P x T-, P-, T- and 1 (no symmetry)-invariant 
solutions, simple classification is impossible, because there are various possibilities of the 
angular momentum dependences of the Q. This situation demands us futher investigation 
in more detail in the future. It must be noted that we have to solve QED-HF equations 
to find which type of the density matrices ( and solutions ) will be realized. The natural 
orbital analysis of a resulted solution is useful to elucidate the type and magnitude of broken 
symmetry. 



The standard approach for inclusion of fluctuation from the HF mean field approximation 
is the MBPT approach [6]. For the purpose, the total Hamiltonian in (9) is devided into 
the zeroth-order and interaction parts as 



where H\ ' denotes the HF mean field interaction potential. The remaining interaction 
Hint is regarded as the perturbation to the HF solution in the case of the MBPT approach [6] . 

However, such single-reference (SR) MBPT approach often breaks down in the case of 
near degenerate ground states, for which generalized QED-HF solutions exist, because of 
the HF instability in (27) [11]. To consider these problems, we introduce the density matrix 
which is defined as 



IV. POST QED-HF METHODS 



A. QED-CASSCF and CASPT2 methods 



Ho = Hk + Hj HF ^ , 
Hm — Hi — Hj HF \ 



(45) 



(46) 
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Pijin) Pij ( K ) 



(47) 



y ^<$(«)|atSt|$(«)> (stofoSJisM) y 

where we write (a, (3, (3*, 7) as (/c). Then we diagonalize the density matrix of the most 
stable QED-HF solution for systems under consideration [14], we obtain the natural orbital 
given by four- component bispinor rj^ and rrii as its occupation number. By using the 
occupation number mj, we can select the active space for treatments of the near degeneracy 
effects [18]. Then the trial CI wavefunction for the QED-MCSCF are easily constructed. 
The QED-MCSCF wavefunction is obtained like the nonrelativistic case: 



l^r 



£cfe^|<t> 7 >. 



(48) 



Here is the configuration state function ( CSF ), Cj is the CI coefficient, and e lSl is 
the orbital rotation operator given in previous section. With introducing the orthogonal 
complement to (48), \^k) = J2jCj K ^\^j), the variational operator of the MC function is 
given as 



T= E T ^{l^)(^o| + |^o)(^x|}, 



(49) 



where, T' 1 " = f and S\ = Si will be satisfied. Then we can write the transformed energy as 



E = (tfo|e~ if e- <5l ife i5l e <f, |tfo) 



(50) 



This expression is, at a glance, same as the nonrelativistic case. Only the definition of 
the orbital rotation operator is different. From (50), we can formally derive the variational 
condition such as the generalized Brillouin theorem in our QED theory. Therefore the multi- 
reference ( MR ) MBPT approach based on QED-MCSCF may be feasible like CASPT2 [19] 
for the nonrelativistic case. 



B. QED-CC and MRCC methods 



The coupled-cluster wave function [20~22] can also be introduced to obtain dynamical 
correlation corrections, by obeying the exponential ansatz: 
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\*cc) =e iS \%). (51) 

Here the operator S is given as 

S = Y,S n = S 1 + S 2 + --- + S N + ---. (52) 

n 

n denotes the excitation level. |$ ) is the QED-HF state. In the relativistic theory, the 
expansion given above can not be truncated at N-excitation level, because of the presence 
of the Dirac sea. Here, the excitation operator S is approximately truncated in the second 
order; S = Si + £>2- Even under this approximation, S2 includes the operators given as 
follows: 

ajajolai, a\aja\b\, djajbhCii, a\ajbkb\, a\Pja\b\, 

a\b]b k a h aJSJSfcSf, kajhai, kafikb}, htyb k b\. (53) 

With considering that the operators a^a, afb>, ba, bb' in exp(iSi) form the Lie algebra, 
there is no other posibility. To write the exponential form, the excitation operators have 
to be bosonic operators. Thus we have to devide the one-particle function space to two; 
occupied and unoccupied spaces. Quantum number i, j, k, I will satisfy some conditions 
like the nonrelativistic case. These operators have to take the form which will give two 
particle excitation with respect to HF state. The operators given above include redundant 
operators, which will not give the excitation with operating the HF state. The Schrodinger 
equation is given as 

He iS \$ ) = Ee i§ \<Z> ). (54) 

Then the projected coupled-cluster equations are given as 

(%\He l§ \%) = E , (55) 
($o\S n He iS \$ ) = E (%\S n e^\%). (56) 

If we replace |$ ) in (51) with |^ ) in (48), we may have the MRCC formulation [14,23,24]: 
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I^mrcc) = e^o> = e^ 1+ ^|* >- (57) 

However, detailed formulations are rather complex because of redandancy. They will be 
discussed elsewhere. A Schema of this work given in Fig. 2. 

V. DISCUSSION 

In this paper, we have developed our method which combines QED with many-body 
techniques. We have given the QED Hamiltonian written by creation-annhilation opetrators, 
and introduced the relativistic Slater determinant in the Thouless form. We have performed 
the group-theoretical classification of the density matrix. After these preparations, we have 
discussed the possibilities of the MCSCF and coupled-cluster method in QED. 

Now we discuss the relation between our theory and other relativistic theories. 

In the context of the electronic structure of atoms, there are three effects in atoms: The 
electron correlation, the relativistic effect and the QED effect [25] . The electronic structure 
of atoms is determined by the relation of these three factors. The electron correlation 
depends on the electron numbers and can be treated by several many-body techniques. The 
relativistic effect becomes large with increasing the atomic number. The relativistic effect in 
neutral or almost neutral heavy elements can be treated satisfactorily by the Dirac-Coulomb- 
Breit no-sea scheme [6], as discussed in introduction. 

In fact, our QED Hamiltonian with neglecting positron states and adopting only the 
density-density interaction taken into account, will derive the Dirac-Coulomb Hamilto- 
nian [26,27] ( as illustrated in Fig. 3 ). Futher, we can have the Dirac-Coulomb-Breit 
Hamiltonian with adding the Breit operator. Moreover, if the small component x is approx- 
imated by using the large component 0, and take into account the fact that x gives small 
contribution, then the Dirac equation can be reduced to the Schrodinger-Pauli equation [28], 
for which 2-component spinor ( GSO ) will correctly allow for relativistic effects up to 
iy/c) 2 . These terms will be calculated by ab initio GSO program package [4,5]. 
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On the other hand, the case to describe the inner core electrons of heavy elements or, 
electrons of highly ionized heavy elements such as lithium-like uranium, the QED effect can 
not be neglected and we must take the Dirac sea into account. This effect can be treated 
by the perturbation theory in QED [9]. But this method can only take into account the 
one-particle QED correction. Practically this method treat the system of atoms with a few 
electrons. For an illminative example of the QED effects, Fig. 4 shows the collision of two 
uranium atoms [29] . This collision generates the pair-creations, and the QED effects might 
be appeared clearly. Moreover, the near orbital degeneracy effects will become important in 
this case, and also in clusters of these atoms or ions. 

In the case of heavy atoms, as the ionicity becomes high and the electron number de- 
crease, the many-body effect becomes small and the QED effect becomes large. Therefore 
we propose that our theory should be applied to the cases where both the many-body effects 
and the QED effects can not be neglected. Heavy elements in middle level of ionicity should 
be one of the subject for our theory. The Dirac-Coulom-Breit no-sea scheme has been ap- 
plied to only neutral or almost neutral atoms. The method of the QED correction works 
only for the highly ionized heavy atoms. Our theory should be suitable for intermediate case 
between them. The collision of two uranium atoms is also one of the interesting subject of 
our theory. Recently, experiments of x-ray irradiation to cluster plasma are performed, and 
verious new phenomena were studied [30] . The cluster of heavy ions under middle revel of 
ionicity is now an interesting subject in this area. The importance of the relativistic and 
QED effects is discussed in such objects. Fig. 5 illustrate scope and limitation of these the- 
ories. It is noteworthy that many interesting phenomena appear in the intermediate regime 
even in the nonrelativistic case [31]. An application of the QED scheme to treat the strong 
electron-electron interaction in solid state physics is also discussed in relation to spin density 
wave and charge density wave states of mixed valence ( MV ) systems [32]. 
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